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A Universal Scaling Law for Jets of Collapsing Bubbles 
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Cavitation bubbles collapsing and rebounding in a pressure gradient Vp form a "micro-jet" en- 
veloped by a "vapor jet" . This letter presents unprecedented observations of the vapor jets formed 
in a uniform gravity- induced Vp, modulated aboard parabolic flights. The data uncovers that 
the normalized jet volume is independent of the liquid density and viscosity and proportional to 
C, = |Vp|i?o/Ap, where Rq the maximal bubble radius and Ap is the driving pressure. A derivation 
inspired by "Kelvin-Blake" considerations confirms this law and reveals its negligible dependence of 
surface tension. We further conjecture that the jet only pierces the bubble boundary if ^ > 4 • 10~^. 

PACS numbers: 47.55.dp,47.55.dd,43.25.Yw 



Jets produced by cavitation bubbles play a key role 
in cutting-edge technologies [l|-|3|] and erosion 
These jets typically arise when a bubble collapses in 
a liquid of anisotropic pressure: At the ending collapse 
stage, the bubble surface develops a fast (>100ms~^ 
[4, 7, 8]) liquid jet. This "micro-jet" is directed inwards 
against the local pressure gradient [9|, defined in 
the absence of the bubble. While the bubble bounces 
off its enclosed gas, the micro-jet pierces the bubble 
and starts penetrating the liquid 0-0, |^ unless hit- 
ting a boundary. During the regrowth ( "rebound" , see 
Fig.fTJi) of the bubble, the micro-jet becomes visible be- 
cause of its conical shell of vapor (Bl, fTll - [T3| , here called 
the "vapor jet" (Fig. [3j). The velocity and structure of 
jets of bubbles were modeled and measured in various 
cases 0, H, m, [13, but no general relation is known 
between the 'jet size' and the underlying pressure gra- 
dient. Such a relation would allow the modulation of 
jets by specifically engineering the pressure field; and, 
vise versa, permit a measurement of a pressure field 
via static images of jetting bubbles (Fig. [1]). 

This letter expands the state-of-the-art in three 
ways: (i) It presents the first high-speed movies of 
the jets caused by a gravity-induced pressure gradi- 
ent ^p = pg in normal gravity {g = 9.81 ms~^, liquid 
density p ~ lO^kgm"^). (ii) It performs a system- 
atic study of the vapor jets observed while varying the 
maximal bubble radius, the liquid viscosity, the liquid 
pressure, and the pressure gradient. The latter is var- 
ied through a modulation of g aboard parabolic flights 
[25| . (iii) A statistical analysis of the data, backed- 
up by a theoretical derivation, reveals that the jet size 
scales with a dimensionless jet-parameter (. 

Our experiment relies on a gravity-induced Vp, 
which exhibits the unique advantage of being uniform 
in space and time. Such a gradient approximates, to 
first order, any smooth field p{x) = p{0) + x^^p + 
0{d'^p)^ where = Vp(0). Examples of are 

pg gravitational field, (la) 

—p{v-^)v stat. potential flow, (lb) 

^0.2RoAph/h^ rigid flat surface, (Ic) 

-0.2RoAph/h^ free flat surface, (Id) 

where v is the velocity field, Rq is the maximal bubble 
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FIG. 1: Observations of the vapor jets directed against 
Vp during the rebound of cavitation bubbles, (a-d) match 
the cases of Eqs. ([T^-d): (a) video on-line Collapse and 
rebound of a bubble {Rq ^ 4 mm, Ap ^ 15kPa) in 'Og' 
(upper) and 'Ig' (lower); note the shock at the collapse, (b) 
Rebounding bubble {Rq ^ 1 mm, Ap ^ 100 kPa) moving 
leftwards while jetting against the dynamic Vp, orthogonal 
to the calculated p-contours. (c, d) Bubbles {Rq ^ 2 mm, 
Ap 100 kPa) rebounding with a jet towards a flat rigid 
surface (h = 5.3Ro) and away from a flat free surface (h = 
5.1i?o), respectively. Images were taken using the setup of 
this letter (a,c,d) and a cavitation tunnel J^] (b). 



radius before collapse, h is the shortest vector from the 
surface to the bubble center, and Ap = po — p^ with 
Po being the pressure at cavity level and Pv the vapor 
pressure. Eq. (|lb[) follows from momentum conserva- 
tion of an incompressible, stationary potential flow. 
Eqs. ([Tel fTd|) provide 'effective' time-averages of the 
self-generated pressure anisotropy of a bubble grow- 
ing and collapsing close to a flat boundary: Vp is 
defined such that adding a constant counter-gradient 
— Vp would suppress the jet according to Ref. [16]. 
Jets described by Eqs. ([T^-d) are shown in Fig. [H 

Experimental setup - Our setup (Fig. [2]) uses a high- 
speed camera (Photron SAl.l) operating at up to 
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250,000 fps with exposure times of 370 ns to record a 
cavitation bubble generated by a pulsed laser (Quan- 
tel CFR 400, 532 nm, 8 ns). The laser is focussed 
inside a liquid volume to form a point-plasma fv^ 
(diameter ^0.1 mm), which quickly cools and conden- 
sates while growing a bubble that subsequently col- 
lapses and rebounds. Whereas past studies [El, [l, [HI 
used lenses to focus a laser, we here use for the first 
time a concave parabolic mirror (Fig.[2j). We found the 
mirror technique to produce bubbles of much higher 
sphericity, since reflection is independent of the liq- 
uid's refractive index and mirrors allow large angles of 
convergence (here 53°) without spherical aberration. 
Our millimeter-sized bubbles are so spherical that the 
tiny gravity-induced pressure difference between their 
top and bottom becomes the dominant source of jet 
formation. To our knowledge, this experiment pro- 
vides the first clean movies of gravity-jets in normal 
gravity conditions. Similar observations in the past 
required large bubbles {Rq > 1 cm) in hyper- gravity. 

The four controllable experimental parameters are 
the maximal bubble radius Rq (varied in the range 
1 — 7 mm), the liquid pressure at cavity level po 
(8 — 80kPa), the norm of the pressure gradient |Vp| 
(0 — 18kPam~^), and the dynamic viscosity 77 (1 — 
30mPas). These parameters are controlled as follows: 
A pressure-regulated vacuum pump depressurizes the 
liquid at a precision of 0.2 kPa, while also remov- 
ing traces of laser-generated gas. The flight manoeu- 
vres (93 ballistic trajectories, straight cruise, 24 steep 
turns) provide intervals of stable gravity at 'Og', 'Ig', 
'1.2g', '1.4g', '1.6g', and '1.8g' (i. e. g = 1.8-9.81 ms-^), 
as well as transition phases, thus offering a wide range 
of gradients | Vp| = pg. By adjusting the energy of the 
laser pulse and the pressure po, bubbles of various radii 
Ro can be obtained. Rq is then measured at 10 /im- 
accuracy on the high-speed movies (e. g. Fig.[T^, left). 
These movies resolve the initial growth and collapse of 
the bubble into more than 100 frames. Demineralized 
water is used in the experiments at variable gravity. 
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FIG. 2: Schematic view of the experiment flown on 
parabolic flights. The test-chamber, filled with a water- 
glycerine mixture of adjustable viscosity, is pressure- 
controlled. An ultra-spherical cavitation bubble is pro- 
duced by a 8 ns laser pulse focused with a parabolic mirror 
at a high convergence angle of 53°. (Dimensions in mm) 




(c) 



1 

d 1 



















(d) 




FIG. 3: (a, b) Observations of the gravity- driven jet of a re- 
bounding cavitation bubble (i?o = 3 mm, Ap = lOkPa) in 
normal gravity: (a) using a back-illumination (see Fig. [2}, 
(b) using a front-illumination and adaptive overlaying of 
different exposures to increase the dynamic range and 
sharpness. The vapor jet envelops a narrow micro-jet in 
agreement with simulations [14]. (c-e) Model (see text). 



while ground-based follow-up experiments use water- 
glycerol mixtures to expand the viscosity range from 
7^ = 1 mPas (pure water) to 77 = 2 mPas (25% glycerol 
mass) and 77 = SOmPas (75% glycerol). The addition 
of glycerol mainly alters 77, but it also affects Pv, and 
the sound speed c. These variations are accounted for 
in the data analysis below (Eq. [3]). 

Experimental results - Observations of rebounding 
bubbles with a gravity-jet are shown in Figs. [1^, [3^, 
[3]3, m Most of these images use a parallel back- 
light (Fig. [2j), hence the bubbles appear in absorp- 
tion against a white background. This type of imag- 
ing allows a precise measurement of the bubble geom- 
etry and a visualization of shocks (Fig. [Iji). An alter- 
native front-illumination discloses the internal struc- 
ture, namely the narrow micro-jet inside the vapor jet 
(Fig. (3)3). The missing jet in 'Og' (Fig. [1^, top) proves 
the gravitational origin of the jets. Ray leigh- Taylor in- 
stabilities during the rebound [19,] can be excluded as 
jet-drivers because of the torus-topology (Fig. [3]3, c). 

Nearby boundaries can also cause jets by altering 
the pressure gradient [sl, [gI, [ii| , and they dominate over 
gravity if A = h'^pg/{RoAp) < 0.2 [square of Eq. 8.8 
in Ref. |l6|. Here, h = 55 mm is the distance from 
the parabolic mirror to the bubble centre. To guaran- 
tee accurate results we only retain gravity-dominated 
cases with A > 0.5, thus keeping a sample of 104 bub- 
bles with jets. Yet, many small bubbles {Rq ^ 2.5 mm, 
thus h/Ro > 22) in 'Og' (A = 0) yield no jet (Fig. [I^, 
top), since the influence of boundaries is too weak. Al- 
beit excluded from the analysis, these data is shown in 
Fig. m as a single point. 

For each bubble in the sample, each high-speed im- 
age of the rebound phase is decomposed into a circu- 
lar disk and a jet, using a x^-^^ of a circular top-hat 
function. The jet volume V-^^ (vapor+liquid. Fig. [3t) 
is then calculated assuming axial symmetry about the 
jet-axis. Since V^l^ only contains a part of the micro- 
jet, we define an effective volume Vjet as the geomet- 
rical extension of V-^^ into the bubble (Fig. [3t) - an 
approach justified in § 'theoretical model'. The rela- 
tion between Vjet and V^l^ depends on the cone angle cp. 
This angle is measured along the edges of the vapor jet 
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rather than at its tip to bypass potential deformations 
of the tip by surface tension. When V-^^ is maximal, 
we observe (p ^ 4° across all bubbles. Trigonometry 
then implies [26| V^^'''' = {V*^''''^^^ + 0,2Rif where 
Ri is the maximal bubble radius during the rebound. 
We finally define a "normalized jet volume" as 

ejet ^ T^r7[(4V3)i??]. (2) 

Aiming for a model of ejet, we adopt the Ansatz that 
ejet is proportional to a non-dimensional parameter 
defined as a power law of the parameters Rq^ p, 
Po, the liquid compressibility (sound speed c). 

The surface tension a is neglected as justified for our 
relatively large bubbles [l^, but we will show (theory 
below) that even the jets of much smaller bubbles re- 
main insignificantly affected by surface tension. The 
most general non-dimensional form of ( then reads 

C = {R^ / g Ap—^-i rj-^+^y , (3) 

where a, /3, 7 are free parameters. To determine a, 
/3, 7 we perform a X^-^^? minimizing the uncertainty- 
weighted rms of Cjet/C the 104 data points. This 
yields a = 1.04 ± 0.03, /3 = 1.05 ± 0.20, 7 = 0.98 ± 0.10 
(with = 0-9), where the ranges are 67% confi- 
dence intervals obtained by bootstrapping the data 
[i^. Since a = /3 = 7 = lis consistent with the 
data, it seems natural to adopt this choice. Substitut- 
ing |Vp| = pg then reduces Eq. (|3]) to 

C = \Vp\Ro/Ap. (4) 

Fig. m shows the measured values of ejet as a function 
of C together with the linear regression (solid line) 

ejet = 5.4 C. (5) 

A remarkable feature of this proportionality relation 
(Eqs. m [5|) is its independence of the viscosity 7^, as 
verified for 77- variations by a factor 30 (Fig. [4j). 

If no vapor jet is observed {V^^^"^ = 0), Eq. ([2]) 
implies ejet = e^^j" ^ 0.002. The inequality e]^j" > 
refiects that when no vapor jet forms, a micro-jet may 
still be present. In fact, no vapor jet arises in all cases 
between no micro-jet (Fig. [3li) and a micro-jet that 
just touches the bubble surface (Fig. [3^). This range 
of indeterminacy is shaded in Fig. HI Using Eq. (j5j) the 
threshold value ("c, where ejet(Cc) = Cj™^ reads ("c ^ 4 • 
10~^. The jet only pierces the bubble surface if > Cc 
i. e. Cc delimits the topological transition between a 
sphere and a torus. To confirm this conjecture further 
investigations of this transition are needed. 

Theoretical model - The relation ejet oc ( will now 
be derived from first principles. By conservation of 
momentum, the micro-jet momentum equals the in- 
tegrated momentum accumulated by the liquid dur- 
ing the bubble growth and collapse. This momen- 
tum, called "Kelvin impulse", was explored by Blake 
[l3|. It can be computed as / = J^^ dt J^^^^dF^ 
where Tc is the collapse time, S is the bubble surface 
(here assumed spherical), and dF = —pdS is the force 
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FIG. 4: Scaling law. Black points are data at varying i?o, 
po, rj and fixed g = 9.81 ms~^. Colored points are data 
at varying Ro, po, g and fixed 77 = ImPas (green: g < 
llms~^, red: ^>16ms~^, blue: intermediate). Some 67% 
measurement uncertainties are shown by the error bars. 
The solid line is the weighted regression (Eq.[5]). The zone 
of experimental indeterminacy, covering the cases between 
Fig. [Sjl and Fig. [3^, is grey-shaded. 



acting on the bubble surface. Spherically symmetric 
(isotropic) terms in the pressure field p vanish in the 
integral over S. Hence, we only consider anisotropic 
pressure terms, here given by a constant Vp, defined as 
the pressure gradient in the absence of a bubble. The 
bubble adds an additional gradient that varies with 
the normalized time r = t/Tc. This bubble-generated 
gradient contains a radial, spherically symmetric term 
(§3.2.3 in [21]) vanishing in the momentum integral, 
and a linear term proportional to caused by the 
motion of the bubble center. Thus, neglecting isotropic 
terms, dF= —f{T)(Vp ■ R) dS^ where /(r) is a scalar 
function and R{t) is the vector from the bubble cen- 
ter to a surface element. In this generic model, the 
micro-jet momentum I solves to ^281] 

I^-VpRlT, . (6) 

Tc is the key term, where side-effects can inter- 
vene. While the Rayleigh-theory [22] implies Tc ~ 
0.915Ro^/p/Ap, the Plesset-theory [23] details that: 

(i) Incondensable gas increases Tc (see § 'Discussion'). 

(ii) Surface tension decreases Tc - by less than 1% 
for our bubbles and by about 7% for a micro-bubble 
{Ro = 10 jam) in water at standard conditions (a = 
0.07 Nm-\ po = 10^ Pa). Carried along to Eq. pp]) . 
this 7%-effect reduces the jet mass by only 14%, com- 
parable to a 1 /im- measurement error of Rq. (iii) Vis- 
cosity increases Tc for small bubbles, but the effect is 
even weaker - about 1% for a bubble with Rq = 10 /im 
(at 77 = 1 mPas, po = 10^ Pa). Hereafter Tc in Eq. ([6j) 
will hence be approximated as 0.915Ro^/p/^p^ 

In analogy to the "Kelvin impulse" we now introduce 
a kinetic "Kelvin energy" T; = 2 /^^^^ \dF - di^|, 
resulting from the work done by the same anisotropic 
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forces that generate the jet's momentum. Thus [29], 

E cx I Vp| Rt . (7) 

Defining m and v as the mass and spatiahy averaged 
velocity of the fuhy developed micro-jet (Fig. [3]) implies 
/= mv and E oc mv^ . Eqs. (j6l[7]) then yield 

V (X — Y^Ap/p e, (8a) 
mcx \Vp\R^p/Ap, (8b) 

where e = Vp/|Vp|. Eq. (|8a|) is a known relation 
[H, while Eq. (|8bp is of interest regarding the jet vol- 
ume. We hypothesize that the effective jet volume 
scales with the micro-jet volume, V^^^^ = emp~^^ at 
an efficiency e oc R\/R^. This approximation derives 
from the observation that the vapor jet grows out of 
the rebound bubble, thus consuming a fraction of the 
rebound volume ex R\. Eqs. ([2| I8bp then imply 

ejet(x|Vp|i?o/Ap = C. (9) 

Alternatively, Eq. ^ also results from assuming ejet 
proportional to E/Eo, where Eq = (47r/3)i?QAp. 

Discussion - In summary, Eqs. (j5l [9j) demonstrate 
experimentally and theoretically that the normalized 
jet volume ejet is proportional to Along this dis- 
covery subtle issues were encountered that are worth 



explaining: (i) Eq. ([9]) is not 'scale invariant' as it de- 
pends on Rq: given Vp, big bubbles yield larger ejet 
than small ones. Formal scale invar iance is nonethe- 
less recovered using an adequately normalized gradient 
V = Ef=iS/(S^i/^o)ei ^ C = |Vp|/Ap. (ii) How 
can a micro-jet survive inside the hot gas during 
the collapse point? This feature might be attributed to 
the life-time of the hot gas (< 1 /is, [24]) being too short 
to evaporate the jet. (iii) We neglected incondensable 
gas inside the bubble, since the measured bubble radii 
R{t) agree within 1% with the Rayleigh-equation [22], 
when neglecting incondensable gas. If bubbles contain 
significant amounts of incondensable gas, they do not 
fully collapse, hence less concentrating their Kelvin- 
impulse, (iv) While Eqs. (|5l[9]) rely on stationary Vp's, 
they also approximate non-stationary situations, if the 
characteristic time-scales are comparable to or above 
Tc. This is illustrated by the data of Figs. [It, d (where 
ejet = 0.2 Rlh-'^) plotted as stars in Fig. H (v) We 
here considered jet parameters < 0.008. Larger 
may produce more complex bubble-jet morphologies, 
requiring a decomposition into spherical harmonics. 
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